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Abstract 

The Pierre Auger Observatory is a hybrid detector for ultra-high energy cosmic 
rays. It combines a surface array to measure secondary particles at ground level 
together with a fluorescence detector to measure the development of air showers 
in the atmosphere above the array. The fluorescence detector comprises 24 large 
telescopes specialized for measuring the nitrogen fluorescence caused by charged 
particles of cosmic ray air showers. In this paper we describe the components of 
the fluorescence detector including its optical system, the design of the camera, 
the electronics, and the systems for relative and absolute calibration. We also 
discuss the operation and the monitoring of the detector. Finally, we evaluate 
the detector performance and precision of shower reconstructions. 

Key words: cosmic rays, fluorescence detector 
PACS: 96.40.-z, 96.40.Pq, 98.70.Sa 
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1. Introduction 

The hybrid detector of the Pierre Auger Observatory [l| consists of 1600 
surface stations - water Cherenkov tanks and their associated electronics - and 
24 air fluorescence telescopes. The Observatory is located outside the city of 
Malargue, Argentina (69° W, 35° S, 1400 m a.s.l.) and the detector layout is 
shown in Fig. [TJ Details of the construction, deployment and maintenance of 
the array of surface detectors are described elsewhere 0]. In this paper we will 
concentrate on details of the fluorescence detector and its performance. 




Figure 1: Status of the Pierre Auger Observatory as of March 2009. Gray dots show the 
positions of surface detector stations, lighter gray shades indicate deployed detectors, while 
dark gray defines empty positions. Light gray segments indicate the fields of view of 24 
fluorescence telescopes which are located in four buildings on the perimeter of the surface 
array. Also shown is a partially completed infill array near the Coihueco station and the 
position of the Central Laser Facility (CLF, indicated by a white square). The description 
of the CLF and also the description of all other atmospheric monitoring instruments of the 
Pierre Auger Observatory is available in 

The detection of ultra-high energy (> 10 18 eV) cosmic rays using nitrogen 
fluorescence emission induced by extensive air showers is a well established 
technique, used previously by the Fly's Eye 0] and HiRes [B| experiments. It is 
used also for the Telescope Array Q project that is currently under construction, 
and it has been proposed for the satellite-based EUSO and OWL projects. 

Charged particles generated during the development of extensive air showers 
excite atmospheric nitrogen molecules, and these molecules then emit fluores- 
cence light in the ~ 300 — 430 nm range. The number of emitted fluorescence 
photons is proportional to the energy deposited in the atmosphere due to 
electromagnetic energy losses by the charged particles. By measuring the rate 
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Figure 2: Schematic layout of the building with six fluorescence telescopes. 



of fluorescence emission as a function of atmospheric slant depth X, an air 
fluorescence detector measures the longitudinal development profile 4§>{X) of 
the air shower. The integral of this profile gives the total energy dissipated 
electromagnctically, which is approximately 90% of the total energy of the 
primary cosmic ray. 

For any waveband, the fluorescence yield is defined as the number of photons 
emitted in that band per unit of energy loss by charged particles. The absolute 
fluorescence yield in air at 293 K and 1013 hPa from the 337 nm fluorescence 
band is 5.05 ± 0.71 photons/McV of energy deposited, as measured in 0- The 
wavelength dependence of the yield has been described e.g. in Since a 
typical cosmic ray shower spans over 10 km in altitude, it is important to stress 
that due to collisional quenching effects the fluorescence yield is also dependent 
on pressure, temperature and humidity of the air. 

The fluorescence detector (FD) comprises four observation sites — Los 
Leones, Los Morados, Loma Amarilla, and Coihueco — located atop small 
elevations on the perimeter of the SD array. Six independent telescopes, each 
with field of view of 30° x 30° in azimuth and elevation, are located in each 
FD site. The telescopes face towards the interior of the array so that the 
combination of the six telescopes provides 180° coverage in azimuth. Figure 
[2] shows the arrangement of the telescopes inside an observation site. 

Figure 3 depicts an individual FD telescope. The telescope is housed in 
a clean climate-controlled building. Nitrogen fluorescence light enters through 
a large UV-passing filter window and a Schmidt optics corrector ring. The 
light is focused by a 10 square meter mirror onto a camera of 440 pixels with 
photomultiplicr light sensors. Light pulses in the pixels are digitized every 100 
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nanoseconds, and a hierarchy of trigger levels culminates in the detection and 
recording of cosmic ray air showers. 




Figure 3: Schematic view of a fluorescence telescope of the Pierre Auger Observatory. 



This paper is organized as follows. In section [2] we describe the components 
of the optical system of an individual telescope, and in section [3] we focus on 
the telescope camera. The electronics of a fluorescence telescope and the data 
acquisition system (DAQ) of an FD station are described in section Q] The 
details of the calibration hardware and methods are given in section [SJ and 
the performance, operation and monitoring of the fluorescence detector arc 
explained in section [6l Finally, in section [7] we concentrate on the basics of 
shower reconstruction using the measured fluorescence signal, and in section [8] 
we summarize. 



2. Optical system 

The basic elements of the optical system in each FD telescope are a filter 
at the entrance window, a circular aperture, a corrector ring, a mirror and a 
camera with photomultipliers. The geometrical layout of the components is 
depicted in Fig. [4] 

The window is an optical filter made of Schott MUG-6 glass [13]. This 
absorbs visible light while transmitting UV photons up to 410 nm wavelength, 
which includes almost all of the nitrogen fluorescence spectrum. Without the 
filter window, the fluorescence signals would be lost in the noise of visible 
photons. 

The aperture, the corrector ring, the mirror, and the PMT camera constitute 
a modified Schmidt camera design that partially corrects spherical aberration 
and eliminates coma aberration. The size of the aperture is optimized to keep 
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the spot siz^ due to spherical aberration within a diameter of 15 mm, i.e. 90% 
of the light from a distant point source located anywhere within the 30° x 30° 
FOV of a camera falls into a circle of this diameter. This corresponds to an 
angular spread of 0.5°. In comparison, the FOV of a single camera pixel is 1.5°. 
The light distribution within the spot is described by the point spread function 
(PSF) shown in Fig. 




Figure 4: Geometrical parameters of the FD telescopes. 

The schematic view of the spot size diagrams over the whole FOV is shown 
in Fig. [6j where the rows correspond to viewing angles 0°, 10°, 15° and 20°. 
The columns corresponds to different displacements of the camera off the focal 
plane, by changing the camera-mirror distance from —5 mm to +5 mm with 
respect to the nominal separation. The central position is located at a distance 
of 1 657 mm from the primary mirror. The asymmetric shape of some spots is 
due to vignetting and camera shadow. This picture also shows the sensitivity 
of the telescope PSF to the precision of the adjustment of the distance between 
mirror and camera. 

2.1. Segmented mirror 

Due to the large area of the primary mirror (~ 13 m 2 ), the mirror 
is segmented to reduce the cost and weight of the optical system. Two 
segmentation configurations are used: first, a tessellation of 36 rectangular 
anodized aluminum mirrors of three different sizes; and second, a structure 
of 60 hexagonal glass mirrors (of four shapes and sizes) with vacuum-deposited 



1 The image of the point source at infinity on the focal surface of the optical system is 
commonly called the "spot" in optics, but it may be better known as a "point spread function". 
The size of the spot characterizes the quality of the optical system. 
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0° incidence 20° incidence 





Figure 5: Top: The simulated 3D distribution of the light intensity for spots on the optical 
axis (left) and close to the camera corner (right; distance from the optical axis = 20°). The 
size of the imaged area is 20 X 20 mm. Bottom: Fraction of encircled energy as a function of 
spot diameter for the spot on the optical axis (upper curve; angular distance from the optical 
axis = 0°; this curve corresponds to the distribution shown top left) and for the spot close 
to the camera corner (lower curve; angular distance from the optical axis = 20°; this curve 
corresponds to the distribution shown top right). 



reflective coatings. In both cases, all mirror segments have a spherical inner 
radius of 3400 mm, allowing possible deviations up to 3420 mm. 

The 12 telescopes at the Los Leones and Los Morados sites use aluminum 
mirrors. The mirror elements were produced from 20 mm-thick cast aluminum 
blocks, and were milled to the required spherical inner radius using a diamond 
milling technique. After the initial milling, the reflective surface was created by 
gluing a 2 mm sheet of AlMgSiOs alloy onto the inner surface of each element. 
The sheets were attached to the aluminum block at elevated temperature and 
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Displacement of the camera off the focal plane (in mm) 




c 



Figure 6: Spot diagrams of the telescope. 




Figure 7: Two different segmentation configurations of the telescope mirror: /eft: 36 
rectangular anodized aluminum mirror segments; right: 60 hexagonal glass mirror segments. 



pressure to achieve high stability. The material used is well-suited for high- 
precision diamond milling, and allows the final mirror surface to achieve a 
roughness below 10 nm. Finally, a 90 nm thick aluminum-oxide layer was applied 
to the surface by chemical anodization to provide additional protection. 

The remaining 12 telescopes in the buildings at Coihucco and Loma Amarilla 
use mirrors with glass segments. The thickness of the glass is 12 mm, and it 
is composed of SIMAX, a borosilicate glass of the PYREX type. The choice 
of ultra-thin, lightweight glass was motivated by the need to maintain the 
optical stability of the segments at different temperatures and for different 
inclinations of the individual segments. SIMAX is suitable for machine grinding 
and polishing, and has good mechanical and thermal stability. 
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The reflective layer on the glass mirror segments is composed of Al layer, with 
thickness of 90 nm, covered by one SiC>2 layer with thickness of 110 nm. The 
aluminum is used for its high reflectivity, and the silica for its high mechanical 
resistance. The reflectivity has been measured at several points of each segment, 
and the average reflectivity at A = 370 nm exceeds 90%. The surface roughness 
has been measured mechanically and optically, and has an RMS of less than 
10 nm. 

In the glass mirrors, the mirror segments are grouped by their curvature 
radii to minimize the radius deviations in each telescope. The curvature and 
the shape of the reflective surface of the segments are measured and controlled 



using standard Ronchi and/or Hartmann tests 

The alignment of the mirror segments is accomplished by directing a laser 
onto each segment. The laser is mounted on the center of curvature of the 
full mirror. The laser beam approximates a point source, and each segment 
is adjusted such that the beam is reflected to a common point (the center of 
the laser aperture). This procedure is used to adjust all the mirror segments 
and match the full mirror to the required spherical shape, with the center of 
curvature aligned with the optical axis of the telescope. 

2.2. Corrector ring 

A novel solution of the optical system with a "corrector ring" was designed 
to keep the advantage of a large aperture of the Schmidt system, and 
simultaneously simplify the production of the element, minimize its weight and 
cost, while maintaining the spot size within the limits of aforementioned design 
specification The aperture area of the telescope with the corrector ring 

is almost doubled with respect to the optical system without any correcting 
element. The analysis of real shower data [l3[ has compared the performances 
of FD optical systems with and without corrector ringj, and has verified that 
the corrector ring enhances the FD aperture by a factor of ~ 2. 

The corrector ring is the circumferential part of the corrector plate of a 
classical Schmidt camera with one planar side and the other with an asphcric 
shape corresponding to a 6^-order polynomial curve. Such a surface is difficult 
to manufacture and therefore some optimizations were adopted to simplify 
the lens production. Eventually, a spherical approximation specially designed 
for the fluorescence detector was chosen to fulfill both price and performance 
requirements (see Fig. [8j. The simplified corrector ring is located at the aperture 
and covers the annulus between radii of 0.85 m and 1.1 m. 

Since the rings have an external diameter of 2.2 m, their manufacture and 
transportation to the site in a single piece would have been very difficult. 
Therefore, each lens was divided into 24 segments. The size and profile of one 
segment is shown in Fig.[9l The production of the segments has been performed 



by Schwantz Ltd lj| after assembling a machine with a circular base to hold 



2 Test measurements without corrector ring were realized on several FD telescopes prior to 
corrector ring installation. 
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t = 12 mm 
dJ = 1700 mm 
d2 = 2200 mm 
R ~ 8380 mm 
y ~ 795 mm 



Figure 8: Corrector ring profile using a spherical approximation, the optical axis of the 
telescope is identified by the dash-dotted line with an arrow on the left and d2 is the diameter 
of the aperture. 



the segments, and a disk with diamond abrasive cylinders for the grinding of the 
glass (BK7 glass from Schott [l(|) to the desired profile. After production of 
each ring, its segments were tested for proper shape. To scan the ring profiles, 
a laser beam was pointed towards the curved and flat surfaces and the positions 
of the reflected light from each surface were measured fl2| . 

2. 3. Simulation of the optical system 

To evaluate the overall optical performance of the detector and to confirm 
the theoretical expectations, a dedicated complete simulation of the optical 
system was developed using Geant4 15, a Monte Carlo toolkit for the 



simulation of radiation and light propagation. The tracking of optical photons 
includes refraction and reflection at medium boundaries, Rayleigh scattering and 
bulk absorption. The optical properties of a medium, such as refractive index, 
absorption length, and reflectivity coefficients, can be expressed as functions of 
the wavelength. The application includes a detailed description of the different 
detector elements - UV filter, corrector lens, mirror and camera. The optical 
properties of all materials, such as the absorption length and the refractive 
index, were implemented as a function of wavelength. 

The simulation confirms that all design specifications of the FD optical 
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Figure 9: Technical drawing of one corrector ring segment. The dimensions are given in 
millimeters. 



system are met, i.e. even in the corners of the camera 90% of the light from a 
distant point source is concentrated within a diameter of 15 mm. 




Figure 10: Ray tracing simulation of the optical system of the telescope made using Geant4. 



3. Telescope camera 



The camera [17[ is the sensitive clement of a telescope. It is composed of a 
matrix of 440 pixels located on the focal surface of the telescope. The camera 
pixels are arranged in a matrix of 22 rows by 20 columns (Fig. [Tib'). The 
corresponding held of view is of 30° in azimuth (full acceptance of one row) and 
28.1° in elevation (full acceptance of one column). 
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A cosmic ray shower is imaged on the camera as a line of activated pixels 
having a track-like geometrical pattern and also a clear time sequence. Each 
pixel is realized by a photomultiplicr with a light collector. 

3.1. Geometry 

The pixel array lies on the focal surface of the optical system, which is a 
sphere of 1.743 m radius. The pixels are hexagonal with a side to side distance of 
45.6 mm, corresponding to an angular size of 1.5°. The pixel centers are placed 
on the spherical surface following the procedure outlined in Fig. 1111 where we 
define AO = 1.5° and A0 = 0.866°. The first center is placed at +A6/2 with 
respect to the telescope axis, which is taken as the z axis in Fig. [TT] The other 
pixel centers are obtained with increasing (or decreasing) AO steps. In this 
manner, a row of twenty pixels (corresponding to 30° in azimuth) is built. 

The following row of pixels is obtained by a rotation of A<j) around the x axis. 
To produce the correct staggering between rows, the pixel centers are moved by 
A9/2 with respect to their positions in the previous row. 



y 




b) c) 

Figure 11: Geometrical construction of the FD camera; a) pixel centers are placed on the 
spherical surface in steps of A9 = 1.5° and Arp = 0.866°, b) positioning of the pixel vertices 
around the pixel center, c) the camera with 440 pixels arranged in a 22x20 matrix. 

Once the pixel centers have been defined, each pixel hexagon is determined 
by positioning six vertices around the center. The angular positions of the 
vertices are obtained by moving in steps of A8/2 and A0/3 with respect to the 
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pixel center, as depicted in Fig.lTTb. Equal steps in angle produce different linear 
dimensions depending on the pixel position on the spherical surface. Therefore, 
the pixels are not exactly regular hexagons, but their size varies over the focal 
surface by 3.5% at most. Differences in the side length are smaller than 1 mm, 
and are taken into account in the design of the light collectors (see section 133)) 
and in the analysis. 

3.2. Mechanics 

The camera body was produced from a single aluminum block by a 
programmable milling machine. It consists of a plate of 60 mm uniform 
thickness and approximately rectangular shape (930 mm horizontal x 860 mm 
vertical), with spherical outer and inner surfaces. The outer radius of curvature 
is 1701 mm, while the inner radius is 1641 mm. Photomultiplicr tubes are 
positioned inside 40 mm diameter holes drilled through the plate on the locations 
of the pixel centers. Small holes in the camera body at the pixel vertices arc 
used to secure the light collectors in place. A picture of the camera body is 
shown in Fig. [T^l 




Figure 12: Picture of the camera body. 



Within a finished telescope, the camera body is held in place by a simple 
and robust two-leg steel support made of 5 cm wide C-shaped steel beams. The 
obscuration of the mirror FOV due to the camera support is less than one-tenth 
of that of the camera. Power and signal cables run inside the C-shaped legs of 
the support without producing additional obscuration. 
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3.2.1. Mechanical precision and alignment 

To align an FD telescope, the pixel surface must be placed at the correct 
longitudinal distance from the center of curvature of the mirror, while the 
camera body should be centered on and perpendicular to the telescope axis, 
with the top and bottom sides of the camera parallel to the ground. Ray 
tracing computations have shown that the spot size increases by about 10 % 
when the longitudinal distance from the center of curvature changes by ±2 mm. 
Therefore, the accuracy on the longitudinal position of the pixels on the focal 
surface should be better than ±2 mm. 

The intrinsic accuracy of the rigid metal frame of of the camera body is 
very good, at the level of ±0.1 mm. The point-to-point internal accuracy of the 
pixels on the camera body is at the level of ±1 mm in both the longitudinal and 
transverse directions, due to the positioning of the photomultipliers and of the 
light concentrators on the rigid frame of the camera body. 

The fluorescence buildings are surveyed by standard topographic methods 
and for each bay two accurately determined reference points are marked on the 
floor defining a line which corresponds to the azimuth of the telescope axis. The 
center of the camera is placed in the vertical plane containing this line and at 
the right nominal height. The camera is then aligned horizontally using a digital 
level-inclinometer with a precision better than 0.05° which is placed on the top 
of the camera body. Once the center of the camera is correctly positioned and 
the camera is horizontal, the right orientation is achieved measuring the distance 
of the four corners of the camera body to the mechanical reference point which 
is located at the mirror center. These measurements are done with a commercial 
laser distance-meter with a precision of ±1 mm mounted on the mirror center. 
The alignment procedure should provide a positioning with space accuracy at 
the level of ±lmm and angular precision at the level of one millirad i.e. better 
that 0.1°. 

Measurements of the image of bright stars on the camera focal surface have 
verified that the alignment procedure for mirror elements and camera body 
meet the design specifications [IH . The precision of the absolute pointing of the 
telescopes has been checked to within an accuracy of 0.1 °. 

3.3. The pixel array 

The hexagonal photomultiplicr tube (PMT), model XP3062 manufactured 
by Photonis [19( is used to instrument the camera. Although their hexagonal 
shape represents the best approximation to the pixel geometry, a significant 
amount of insensitive area is nevertheless present between the photocathodes. 
In fact, some space between the PMTs is needed for safe mechanical packaging 
on the focal surface; moreover, the effective cathode area is smaller than the 
area delimited by the PMT glass envelope. To maximize light collection and 
guarantee a sharp transition between adjacent pixels, the hexagonal PMTs are 
complemented by light collectors. 

The basic element for the pixel light collector is a Mercedes star, with three 
arms oriented 120° apart, which is positioned on each pixel vertex. Six Mercedes 
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stars collect the light for a given pixel. The geometrical structure of the light 
collector for one pixel is shown in Figll3h.. 



E i 2 y = SO mm 




x (mm) 

a) b) 

Figure 13: a) Six Mercedes stars positioned in order to form a pixel. Each Mercedes star has 
three arms oriented 120° apart. In the drawing the bottom Mercedes is slightly displaced 
for clarity, b) Measurement of the light collection efficiency, with a light spot moved along 
a line passing over three adjacent pixels. The full (open) dots represent the measurements 
performed with (without) Mercedes stars. 

The Mercedes stars are made of plastic material covered by aluminized 
Mylar™ foils. The arm length is approximately half of the pixel side length. 
The arm section is an equilateral triangle. The base length of 9.2 mm is designed 
to cover the insensitive space due to photocathodc inefficiency (ss2 mm for each 
adjacent PMT) plus the maximum space between the glass sides of the PMTs 
(of the order of 5 mm). The triangle height is 18 mm, and the corresponding 
angle at the vertex is 14.3°. 

Each Mercedes star is held by a bar, about 10 cm long, which is inserted into 
a 3.2 mm hole located on the position of the pixel vertices and kept in place by 
a small O-ring. The geometry of the light collectors was designed on the basis of 
the optical system properties. The range of angles of incidence for the rays on 
the camera is in the interval between approximately 10° to 35°, as determined 
by the shadow of the camera and the aperture of the diaphragm. Note that the 
pixels are defined on the focal surface at 1743 mm from the center of curvature 
where the top edges of the Mercedes divide one pixel from its neighbors. The 
PMT cathodes are recessed behind that focal surface by 18 mm. 

Monte Carlo ray tracing has shown that the light collection efficiency, 
averaged over the camera, is 94%, assuming a reflectivity of 85% for the 
aluminized Mylar foils of the light collector surface. Without the light collectors, 
the collection efficiency decreases to 70%. 

The Monte Carlo simulation was experimentally checked using a specially 
designed light diffusing cylinder with an exit hole for light rays having the same 
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size as the spot produced by the telescope optics. The distribution of the angles 
of the light rays from this optical device is similar to that produced by the 
telescopes. Measurements were made with and without the light collectors. 
Results of a scan moving the light spot over the camera are shown in Fig. [T3"b . 
Without light collectors, a significant loss of light at the pixel borders was 
observed. When installed, the light collectors efficiently recuperate the light 
loss. From these measurements, the light collection efficiency averaged over the 
camera focal surface was found to be 93%. 




Figure 14: a) Detail of the camera body with four PMTs mounted together with two Mercedes 
stars. The large holes to insert the PMTs and the small ones to mount the Mercedes are visible, 
b) Picture of a camera completely assembled with all PMTs and light collectors in place. 

A picture of a small section of the camera and a full picture of the camera 
completely assembled are shown in Fig. [14] 

3.4- The photomultiplier unit 

The XP3062 photomultiplier [l^ is an 8-stage PMT with a hexagonal 
window (40 mm side to side). It is manufactured with a standard bialkaline 
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photocathode with quantum efficiency of about 25 % in the wavelength range 
350-400 nm. The nominal gain for standard operation of the FD is set to 5 x 10 4 . 

The PMT high voltage (HV) is provided by a HV-dividcr chain which forms 
a single physical unit together with the signal driver circuitry. This unit, 
called head electronics (HE), consists of three coaxially interconnected printed 
circuits boards (PCBs): the bias PCB (innermost one), a laser-trimmed hybrid 
technology driver circuit (intermediate), and the interface PCB (outermost). 
The innermost and outermost PCBs are two-sided and of circular shape (32 mm 
diameter), and are interconnected using high reliability pin connectors. The 
HE units were manufactured by Intratec-Elbau (Berlin, Germany) [22| and are 
soldered to the flying leads of the PMT. To ensure central mounting of the HE 
with respect to the symmetry axis of the PMT, and to improve mechanical 
rigidity, a specially-designed plastic structure has been introduced between the 
glass tube and a central guidance hole left in the innermost PCB of the HE. A 
PMT with the attached head electronics is shown in Fig. [TBI 




Figure 15: Picture of a PMT unit. 

The HV divider keeps the PMT photocathode at ground and supplies 
positive HV to the anode. To dissipate minimum power in the dense package of 
HE units, the bleeder is operated at currents of less than 170 /j,A. Stabilization 
of the HV potential for large pulses and in presence of a relatively strong light 
background is thus realized by employing an active network [2fJ that uses 
bipolar transistors in the last three stages of the PMT. This is necessary for 
a fluorescence telescope, which is exposed to the dark sky background and in 
some cases also to a fraction of moon light. The normal dark sky background 
induces an anode current of about 0.8 /iA. The active divider ensures that the 
gain shift due to the divider chain is less than 1% for anode currents up to about 
10 fxA. 

The driver located on the intermediate ceramics board of the HE receives 
the AC-coupled anode signal through a differential input. The anode pulse flows 
through a load resistor and reaches one leg of the differential input of the line 
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driver, while the other leg only picks up the common-mode noise. A common- 
mode rejection ratio of 28 dB is obtained between 1 — 100 kHz. The integrated 
laser-trimmed hybrid circuit supplies a balanced output matching the 120 Q 
characteristic impedance of the twisted pair cable connecting the HE to the 
front-end analog board of the readout electronics (see section l4.1.ip . 



Extensive qualification and acceptance tests [21[ were performed on the full 
PMT unit, i.e. the PMT with the HE attached, using automatic test systems. 
Among the tests were measurements of the HV dependence of the gain, which 
are needed for a proper grouping of the PMTs with similar gain. The relative 
photocathode sensitivity at different wavelengths, the linearity of the PMT 
unit, and the photocathode uniformity were also measured. Before and after 
installation, the relevant information for each individual PMT unit is collected 
into a database which tracks the PMT characteristics, including calibration 
measurements performed in situ. 

To reduce the cost of the power supplies, the photomultipliers of each 
camera are organized into ten groups of 44 units. Each group has similar gain 
characteristics, and is powered by a single HV channel. The spread of the gains 
within a group is about ± 10%. 

Cables are distributed to the PMT electronics through distribution boards 
positioned behind the camera and within its shadow, i.e., without causing 
additional obscuration. These boards serve groups of 44 PMTs of the camera, 
supplying HV and LV and receiving the differential signals from the drivers in 
twisted pair cables. On the board a fuse for overcurrent protection of each LV 
line is provided. From the boards, round shielded cables carry the PMT signals 
on twisted pair wires to the front-end crate, located at the base of the camera 
support. 

4. Electronics and data acquisition system 

The FD telescopes are used to record fluorescence signals of widely varying 
intensity atop a sizeable, and continuously changing, light background. This 
presents a significant challenge for the design of the electronics and data 
acquisition system (DAQ), which must provide a large dynamic range and strong 
background rejection, while accepting any physically plausible air shower. The 
DAQ must also allow for the robust, low-cost, remote operation of the FD 
telescopes. And finally, the absolute FD-SD timing offset must be sufficiently 
accurate to enable reliable hybrid reconstruction. 

The FD electronics are responsible for anti-alias filtering, digitizing, and 
storing signals from the PMTs. As the PMT data are processed, they are 
passed through a flexible three-stage trigger system implemented in firmware 
and software. The remaining high-quality shower candidates are packaged by 
an event builder for offline shower reconstruction. For each shower candidate, a 
hybrid trigger is generated for the surface detector. An overview of the complete 
trigger sequence is presented in Table [T] 

The organization of the electronics and DAQ is hierarchical, reflecting 
the physical layout of the FD buildings. Figure \W\ shows the readout 
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scheme of one FD site, divided into four logical units: head electronics for 
440 PMTs x 6 telescopes, which provide low and high voltage; front-end (FE) 
sub-racks, where the signals are shaped and digitized, and where threshold 
and geometry triggers are generated; the DAQ subnet, in which six telescope 
PCs (MirrorPCs) read out the stored data and perform additional software- 
based background rejection; and an FD site network, in which a single EyePC 
merges triggers from the six telescopes and transfers them to the Observatory 
Central Data Acquisition System (CDAS) in Malargiie. The FD site network 
also contains a Slow Control PC to allow for remote operation of the building. 
The camera electronics and PCs are synchronized by a clock module based on 
the Motorola Oncore UT+ GPS receiver, the same receiver used in the SD array. 

4-1- Front-end electronics 

Each FD camera is read out by one front-end sub-rack and an associated 
MirrorPC. The front-end electronics contain 20 Analog Boards (ABs), and each 
AB receives data from a column of 22 PMT channels. The boards are designed 
to handle the large dynamic range required for air fluorescence measurements; 
at the energies of interest for the Pierre Auger Observatory, this means a range 
of 15 bits and 100 ns timing (23|. The sub-racks also contain dedicated boards 
for hardware triggers: 20 First Level Trigger (FLT) boards for pixel triggers, 
and one Second Level Trigger (SLT) board for track identification within the 
camera image. 
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FD Trigger Sequence 


Trigger Level 


Location 


Purpose 


Event Rate 


1 

(FLT) 


FE sub-racks 
(FLT boards) 


pixel threshold 
trigger 


100 Hz pixel" 1 


2 

(SLT) 


FE sub-racks 
(SLT board) 


track shape 
identification 


0.1 - 10 Hz telescope" 1 


3 

(TLT) 


MirrorPCs 
(software) 


lightning 
rejection 


0.01 Hz telescope -1 


T3 


EyePC 
(software) 


event builder, 
hybrid trigger 


0.02 Hz building^ 1 



Table 1: Trigger sequence for FD events. At each telescope, events are selected based on 
channel thresholds (FLT), track shape (SLT), and lightning rejection (TLT). Events passing 
the TLT are merged by an event builder on the FD EyePC. If the event passes further quality 
cuts for a simple reconstruction, a hybrid trigger (T3) is sent to CDAS. 



The Analog Board 

The purpose of the AB is to receive inputs from the head electronics on the 
PMT camera and adapt them for digitization by the analog-to-digital converters 
(ADCs) located on a corresponding FLT board. As shown in Fig.[T71 the analog 
and FLT boards arc physically connected by three 50-pin SMC connectors. The 
combined front-end module measures 367 mm x 220 mm, and is housed in a 9U 
standard crate next to each FD telescope. Every crate contains 20 modules in 
total. 

One AB channel comprises the following logical blocks: receiver, gain stage, 
anti-aliasing filter, and dynamic range adapter. The channel receives input from 
the HE on the PMT camera via a distribution board. Individual pixel enabling 
is performed by a fast analog switch, which is also used to generate an internal 
test pulse upon request from the FLT logic. In each channel, the AB: 

• performs a differential to single-ended conversion of the signal; 

• adjusts the channel gain; 

• applies an anti-aliasing filter before signal sampling; 

• adapts the 15-bit dynamic range to the 12-bit ADCs; 

• provides an injection point for test pulses. 

The telescope PMTs are coarsely gain-matched during the installation of the 
PMT camera to avoid expensive individual HV control. Additional fine-tuning 
of the channel gains is achieved by means of digital potentiometers connected 
in series with a resistor on the feedback loop of each non-inverting gain stage. 
The devices can change individual channel gains up to a factor of about 1.9, 
and allow for gain matching of the channels in the camera to within 0.6%. 
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Figure 17: Photo of the FLT (left) and Analog Board (right): Both boards are connected by 
three 50-pin SMC connectors. A stiffener bar and common front panel provide mechanical 
rigidity. 



Prior to sampling, the PMT signal is processed by an anti-aliasing filter to 
match the 10 MHz digitization rate. A fourth-order Bessel filter with a cutoff 
frequency of 3.4 MHz has been implemented in the AB as a compromise between 
reconstruction error and circuit complexity. The Bessel filter scheme, featuring 
a linear dependence of the transfer function phase with frequency, was selected 
after a detailed study of optimal filters to minimize distortion of the current 
signal shape. 

The final component of an AB channel is the dynamic range adapter. The 
FD records signal sizes between 3 and 10 5 photoelectrons per 100 ns, covering a 
dynamic range of 15 bits. Rather than digitize the full range of signal sizes, an 
optimal cost-effective solution using dynamic range compression to 12-bit ADCs 
has been implemented in the analog electronics design. 

The compression technique, which uses virtual channels, relies on the fact 
that the the shower signal does not appear on all pixels simultaneously; instead, 
each pixel is triggered in a well-ordered time sequence. Every channel is 
configured with a high and low gain of about 20 and 1, respectively. The 
high gain is optimized for the most frequent small- and intermediate-amplitude 
pulses, and is digitized pixel by pixel. The signal of the low gain stage is routed 
together with signals from 10 other non adjacent pixels in an analog summing 
stage. The sum signal is then processed by a virtual channel with a gain near 
unity. 

When a nearby high-energy shower is detected by the FD, typically only one 
channel out of the group of 11 will saturate its 12-bit range at a given time. 
Even in case of multiple saturation in the same group, this will never occur in 
overlapping time bins. Therefore, the signal can be recovered from the virtual 
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Figure 18: Simulated air shower of 10 21 eV at 10 km from the telescope. The top-left figure 
shows the light track in the cameras. The lower left figure shows the reconstructed shower 
profile without use of the virtual channel information. Making use of the virtual channel (lower 
right), the shower parameters can be reconstructed with minimal systematic distortions. 



channel without ambiguity. This is demonstrated in Fig. 1181 which depicts a 
simulated 10 21 eV shower with a landing point 10 km from the FD telescopes. 
The high-gain channels saturate during the development of the shower, and an 
attempt to reconstruct the shower using only these channels leads to significant 
systematic distortions in the shower profile. However, with the information 
in the virtual channel, the shower is reconstructed with almost no systematic 
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4-1-2. First Level Trigger (FLT) Module 

The First Level Trigger (FLT) module is the heart of the digital front-end 
electronics. The module processes the data from one 22-channel column. Its 
main tasks are to: 

• continuously digitize signals from the AB at 10 MHz; 

• store the digitized raw data in memory for later readout; 

• measure the pixel trigger rate for each channel; 

• compensate for changing background conditions and maintain a pixel 
trigger rate of 100 Hz by adjusting the pixel trigger thresholds; 
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• allow access to raw data memory and internal registers; 



• provide a digital interface to the AB to generate test pulses and to set the 
analog gain at the AB; 

• calculate the baseline offset and its fluctuation averaged over a 6.5 ms 
period; 

• calculate the multiplicity (number of triggered pixels) in one column. 

The functions of the FLT are implemented in FPGA (Field Programmable 
Gate Array) firmware to improve the flexibility and cost-effectiveness of the 
module. A pipelined 12-bit ADC (ADS804) is used to digitize the signal at 
10 MHz in each channel, and the data are stored with a 4-bit status word in 
64k x 16-bit SRAMs. The address space of each SRAM is divided into 64 pages 
of 1000 words. In the absence of triggers from the SLT module (described in the 
next section) , each page works as a circular buffer to hold the ADC values of the 
previous 100 /is. When an SLT trigger occurs, all FLT boards synchronously 
switch to the next unused memory page, whose address is provided by the SLT. 
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Figure 19: Generation of the First Level Trigger in one channel. A moving sum of 5 to 16 
values is compared to an adjustable threshold. The use of the sum improves the S/N-ratio, 
and allows for regulation of the trigger rate. 



The main task of the FPGA logic is to generate the pixel trigger (FLT) using 
a threshold cut on the integrated ADC signal. The FLT is shown schematically 
in Fig. Q3B A moving "boxcar" sum of the last n ADC samples is compared to a 
14-bit threshold. Here n is a fixed number of time bins which can be chosen in 
the range 5 < n < 16. The threshold is dynamically adjusted to maintain a pixel 
trigger rate of 100 Hz. When the sum exceeds the threshold, a pixel trigger is 
generated. The use of the sum substantially increases the signal-to-noise (S/N) 
ratio for each channel; for n = 16, S/N improves by a factor of 2.8 [26|. When 
the moving sum drops below the threshold, a retriggerable mono-flop extends 
the pixel trigger for an adjustable period of 5 /is to 30 /j,s common to all pixels, 
increasing the chance of coincident pixel triggers. 
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The rate of pixel triggers, called the "hit rate," is measured in parallel 
by counters for each channel. It is used to adjust the threshold in such a 
manner that the hit rate is kept constant at 100 Hz under variable background 
light conditions. The background light levels seen by each PMT can also be 
monitored by analysis of the variances of the ADC values. Therefore, the 
FPGAs regularly evaluate the mean and variance of the channel baselines every 
6.5 ms using 65536 consecutive ADC samples (27| |. The measured channel hit 
rates, thresholds, offsets, and variances are periodically queried and stored in a 
monitoring database for experimental control (see section f6.5p . 

Finally, the multiplicity, or the number of pixels triggered simultaneously 
within 100 ns, is calculated for each FLT and for the full camera. The number 
of channels above threshold in one 22-pixel column is tracked on the FLT board. 
The total sum for the full camera is evaluated at the SLT board using daisy- 
chained lines on the backplane. The chronological sequence of multiplicity values 
carries information about the temporal development of the camera image, which 
is used by the software trigger (see section |4 . 2. lj) . 

4-1.3. The Second Level Trigger (SLT) 

The pixel triggers generated for each channel in the 20 FLT boards of the 
FE sub-rack are read out by a Second Level Trigger (SLT) board. The functions 
of the SLT are implemented in FPGA logic, and its primary task is to generate 
an internal trigger if the pattern of triggered pixels follows a straight track [28| . 
The algorithm searches for track segments at least five pixels in length. It uses 
the fundamental patterns shown in Fig. [20l as well as those created by rotations 
and mirror reflections of these segments. 
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Figure 20: Fundamental types of pattern regarded as straight track segments. 



During data acquisition, some tracks will not pass through every pixel center, 
and therefore some PMTs along the track may not record enough light to trigger. 
To allow for this situation, and to be fault-tolerant against defective PMTs, the 
algorithm requires only four triggered pixels out of five. Counting all different 
four-fold patterns originating from the five-pixel track segments in Fig. [201 one 
finds 108 different pattern classes. 

A brute-force check of the full 22 x 20 camera pixel matrix for the 108 pattern 
classes would result in 37163 combinations in total. Therefore, the SLT uses a 
pipelined mechanism that scans for tracks on a smaller 22 x 5 sub-matrix of the 
camera. Every 100 ns, the FLT pixel triggers of two adjacent columns (22 x 2 
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pixels) are read into a pipelined memory structure. During the scan of the full 
camera, coincidence logic identifies the pattern classes 1 to 108 and generates 
a trigger if a pattern is found. The full scan lasts 1 /its, and the pixel trigger 
information scanned during this time is stored in the pixel memory of the SLT, 
which is organized as a ring buffer for 64 events (as on the FLT). The FLT 
multiplicities are also stored in the same memory, providing complementary 
measures of the spatial and temporal size of the event. 

Note that in addition to generating the geometrical trigger, the SLT board 
also supplies bookkeeping information for the FE sub-rack. In particular, the 
board provides a time stamp for each event using synchronization signals sent 
by the GPS module over a dedicated line. 

4- 1-4- Test Systems 

The large number of channels in each telescope and the complexity of the 
electronics and triggers require the use of an automatic testing system. The 
testing and calibration of the full optics-PMTs-electronics chain is described in 
section [5] However, the system also allows for tests of the electronics alone. A 
test-pattern generator has been implemented on the front-end boards to create 
pulses of variable width and amplitude at the input of each analog channel. 
These test pulses are used to check the analog and digital functions of the FLT 
and SLT boards. The system is independent of the full camera, and can be used 
for maintenance when the HV/LV supply of the camera is turned off. 

4-2. DAQ Software and Software Trigger 

Once an event has been processed by the DAQ hardware and stored in FLT 
and SLT memory, it can be read out and analyzed by trigger software in the 
MirrorPCs. Each MirrorPC is connected to the digital boards in the FE sub- 
rack via a Fire Wire interface [24[ , and the MirrorPCs communicate with the 
EyePC through a 10/100 Mbit LAN switch. 

The DAQ system can handle data from three types of sources: 

1. external triggers; 

2. calibration events and test pulses; 

3. real air shower events. 

External triggers are primarily artificial light sources used for atmospheric 
monitoring, such as laser shots from the Central Laser Facility 0| ■ The various 
types of calibration pulses are described in section [5] 

Events triggered externally, or by the calibration and test-pulser systems, 
receive special trigger bits in the SLT status word to distinguish them from 
real air shower data. When the MirrorPCs read out such events from the FE 
modules, they are sent directly to the EyePC for event building and storage in 
dedicated raw data files. In contrast, when the MirrorPCs read out true air 
shower data, the events are processed by a Third Level Trigger. The surviving 
data are sent to the EyePC, which builds an event from the coincident data in 
all telescopes and generates a hybrid trigger (T3) for the surface array. 
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4.2.1. The Third Level Trigger (TLT) 

The Third Level Trigger (TLT) is a software algorithm designed to clean 
the air shower data stream of noise events that survive the low-level hardware 
triggers. It is optimized for the fast rejection of triggers caused by lightning, 
triggers caused by muon impacts on the camera (see Fig. I31|) . and randomly 
triggered pixels. 

During the austral summer, distant lightning can significantly disrupt the 
normal data acquisition of the FD telescopes. In good atmospheric conditions, 
the SLT will detect one to two events per minute per telescope. However, 
lightning can cause large parts of the camera (i.e., hundreds of pixels) to trigger 
in bursts of several tens of events per second. If not filtered, the bursts will 
congest the 64-event buffers in the FLT and SLT boards, and the telescope 
readout (including ADC traces) will suffer a substantial increase in dead time. 

The TLT is designed to efficiently filter lightning events without performing 
costly readouts of the complete ADC traces. The algorithm achieves high speed 
by reading out only the FLT multiplicities and the total number of triggered 
pixels. Using several thousand true showers and background events from one 
year of data, a dedicated study was performed to determine efficient lightning 
cuts based on the time development of the multiplicity and the integral of the 
multiplicity over the whole event consisting of 1000 subsequent multiplicity 
values With these cuts, approximately 99% of all lightning events are 

rejected by the TLT during a short (50 /us) decision time. 

The multiplicity-based lightning rejection effectively removes noise events 
with more than 25 pixels. The remaining events with smaller number of 
triggered pixels, which are contaminated by random pixel triggers and muon 
impacts, are filtered in a second step. The correlation between the spatial 
arrangement and peak signal times of triggered pixels is used to discard noisy 
channels far off the light track [29|. The full ADC traces must be used to 
determine the pixel trigger times, but this readout does not appreciably increase 
the telescope dead time when the number of pixels is small. 

The TLT performance has been validated with simulated showers and data 
recorded during different weather and sky background conditions. Although the 
exact behavior of the algorithm depends on the actual conditions, approximately 
94% of all background events are correctly identified by the TLT. The fraction 
of true showers rejected by the trigger is below 0.7%. 

4.2.2. Hybrid Trigger (T3) 

Events passing the TLT in each telescope are sent to the EyePC through the 
DAQ subnet. An event builder running on the EyePC merges coincident events 
from adjacent telescopes. The EyePC also sends a hybrid trigger, called a T3, 
to the CDAS in Malarguc. 

The T3 acts as an external trigger for the surface array. Its purpose is to 
record hybrid events at low energies (below 3 x 10 18 eV) where the array is 
not fully efficient and would not often generate an independent trigger. Hybrid 
events at these energies occur within 20 km of the FD buildings and usually do 
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not trigger more than one or two SD stations. However, as discussed in section[71 
even limited SD data are sufficient for high-quality hybrid reconstruction. 

The T3 algorithm is used to calculate a preliminary shower direction and 
ground impact time with a simple online reconstruction. Once these data arrive 
at the CDAS, a request is sent to the SD for signals recorded close to the 
calculated impact time. For each T3 trigger, the SD stations nearest the FD 
building (comprising approximately one-quarter of the array) are read out. The 
FD and SD data are merged offline for subsequent hybrid analysis. 

5. Calibration 

5. 1 . Introduction 

The reconstruction of air shower longitudinal profiles and the ability to 
determine the total energy of a reconstructed shower depend on being able 
to convert ADC counts to a light flux at the telescope aperture for each channel 
that receives a portion of the signal from a shower. It is therefore necessary to 
have some method for evaluating the response of each pixel to a given flux of 
incident photons from the solid angle covered by that pixel, including the effects 
of aperture projection, optical filter transmittance, reflection at optical surfaces, 
mirror reflectivity, pixel light collection efficiency and area, cathode quantum 
efficiency, PMT gain, pre-amp and amplifier gains, and digital conversion. While 
this response could be assembled from independently measured quantities for 
each of these effects, an alternative method in which the cumulative effect is 
measured in a single end-to-end calibration is employed here. 

The absolute calibration of the fluorescence detectors uses a calibrated 
2.5 m diameter light source (known as the "drum") at the telescope aperture, 
providing same flux of light to each pixel. The known flux from the light source 
and the response of the acquisition system give the required calibration for each 
pixel. In the lab, light source uniformity is studied using CCD images and the 
intensity is measured relative to NIST-calibrated photodiodes [13] . Use of the 
drum for gain adjustment and calibration provides a known, uniform response 
for each pixel in each camera of the FD detector. The average response of the 
FD is approximately 5 photons/ ADC bin. 

Three additional calibration tools are used as well: 

• Before and after each night of data taking a relative calibration of the 
PMTs is performed. This relative calibration is used to track both short 
and long term changes in detector response. 

• The relative FD response has been measured at wavelengths of 320, 337, 
355, 380 and 405 nm, defining a spectral response curve which has been 
normalized to the absolute calibration. 

• An independent check of the calibration in some phototubes is performed 
using vertical shots from a portable laser in the field. 

The sections below describe the hardware and use of these calibration 
systems. 
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5.2. Light source and drum 

The technique [30l | is based on a portable 2.5 m diameter, 1.4 m deep, drum- 
shaped light source which mounts on the exterior of the FD apertures (see 
Fig. I21[) . The source provides a pulsed photon flux of known intensity and 
uniformity across the aperture, and simultaneously triggers all the pixels in the 
camera. The surface of the drum is a good Lambertian. This means that the 
light emitted per unit solid angle from any small area A depends only on the 
angle 9 with respect to the normal direction according to 1(9) = IqAcos(9). 
Looking at the disk from angle 9 and distance d with a fixed solid angle lu, the 
viewed disk area is dco/ cos(f?). So the intensity is 1(9) = Iodui. The intensity is 
independent of viewing angle for a perfect Lambertian. Fig. [^Ushows that pixels 
of a CCD camera (each seeing a small portion of the drum surface) measure the 
same intensity regardless of the viewing angle to the drum. 

To produce diffuse light inside the drum, illumination is provided by a pulsed 
UV LED [3l| (375 ±12 nm), mounted against the face of a 2.5 cm diameter 
x 2.5 cm long Teflon™ cylinder. The Teflon cylinder is mounted in a 15 cm 
diameter reflector cup, which is mounted flush to the center of the drum front 
surface, illuminating the interior and the back surface of the drum. The LED 
is inserted down the axis of the drum from the back through a pipe. A silicon 
detector attached to the opposite end of the teflon cylinder monitors the light 
output for each pulse of the LED. 

The drum was constructed in sections, using laminations of honeycomb core 
and aluminum sheet. The sides and back surfaces of the drum interior are lined 
with Tyvek™, a material diffusively reflective in the UV. The reflecting surfaces 
of the cup are also lined with Tyvek. The front face of the drum is a 0.38 mm 




Figure 21: A schematic showing the drum mounted in a telescope aperture. 
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thick Teflon sheet, which transmits light diffusively. 

5. 3. Calibration of the drum 

The absolute calibration of the drum light source intensity is based on UV- 
enhanccd silicon photodctectors, calibrated at NIST to ±1.5%. While the small 
surface area and low response of these detectors preclude detection of the small 
photon flux from the drum surface directly, the photodiodc calibration can be 
transferred to a more sensitive PMT/DAQ system. 

To establish the absolute flux of photons emitted from the drum surface, a 
reference PMT is placed on the drum axis, 14 m from the surface. The LED 
light source in the drum is pulsed for a series of 5 fis pulses and the charge from 
the PMT is integrated and recorded for each pulse, resulting in a histogram of 
the distribution of the observed integrated flux. 

On an optical bench, the PMT is then exposed to a small diffuse LED light 
source with neutral density filter. This filter adjusts the light intensity to the 
level similar to that of the drum. The LEDs are pulsed with the same driving 
circuitry as for the drum, and the intensity is set to a series of values producing 
a series of histograms with PMT ccntroids surrounding that from the drum 
measurement above. At each of these intensity settings a second measurement 
is made in which the PMT is replaced by the NIST-calibrated photodiode and 
a neutral density filter in the source is removed, increasing the intensity to a 
level measurable by the photodiode. For this second measurement, the LEDs 
arc run in DC mode. The relationship of PMT response to photodiode current 
is found to be very linear. The flux of photons at 14 m from the drum surface 
is calculable from the active area of the photodiode, the neutral density filter 
reduction factor, the LED pulsed/DC duty factor, the NIST calibration for the 
photodiode, and the value of photodiode current corresponding to the PMT- 
drum centroid. The value of photodiode current is interpolated from the linear 
response-current relationship given above. The resulting uncertainty in the 
drum intensity is 6%. Additional uncertainty contributions related to use of the 
drum at the FDs [35j |. such as temperature dependencies, along with camera 
response variations in time and spectral characteristics of the LED light source, 
combine with the drum intensity uncertainty resulting in an overall uncertainty 
of 9% for the absolute FD calibration. 

5.4- Drum relative uniformity measurements 

Uniformity of light emission from the drum surface is important, since the 
pixels in a FD camera view the aperture at varying angles. In addition, for each 
pixel, a different part of the aperture is blocked by the camera itself. Studies 
were made of uniformity of emission across the surface and as a function of 
viewing angle. These uniformity measurements were made using a CCD, viewing 
the emitting surface of the drum from a distance of 14 m (see Fig. |2"2"| . Images 
were recorded with the drum axis at angles of 0, 10, 20 and 25 degrees relative to 
the CCD axis, covering the range of the Auger telescope field of view (0 — 21°). 
For these images, the UV LEDs were powered continuously. 
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A 0° image was used to analyze the uniformity of emission over the drum 
surface. Using software, concentric circles were drawn on the image, defining 
annular regions of increasing radius, as shown in the figure. The intensity of 
the pixels in each region was analyzed to obtain the intensity as a function of 
radius. In the area defined by the 2.2 m aperture radius, the relative uniformity 
of intensity is constant over the area to about ±2% except for a central dark 
spot. The variation with viewing angle of a section of the drum image is also 
shown in Fig. fJH 
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Figure 22: Left: CCD image at 0° drum angle, showing the defined rings for relative intensity 
analysis. Production deformations in the teflon material can be seen. Middle: A plot of the 
observed pixel intensities in the defined regions of the drum, shown in the left panel; Right: 
The results of angular relative intensity measurements at 0, 20, and 25°. inclinations of the 
drum with respect to the CCD camera. 

The measured drum non-uniformities are small (diagonal stripes in Fig. |2"21 
intensity decrease with radius, etc.), indicating that the FD pixels see similar 
intensities integrated over the drum surface. While perfect drum uniformity is 
desirable, the present non-uniformities are acceptably small and well mapped 
over the surface of the drum. A ray-tracing program using the uniformity and 
angular intensity information from the CCD images shows less than 1% variation 
in total flux seen by the pixels, and corrections are applied for these variations. 



5.5. Absolute calibration using vertical laser shots 

The drum technique for absolute calibration has been checked for some pixels 
using remote laser shots at 337 and 355 nm 34j, [35j. A laser pulse is shot 
vertically into the air with a known intensity. A calculable fraction of photons 
is scattered to the aperture of the FD detector. This yields a known number 
of photons arriving to the detector for each pixel, which views a segment of 
the laser beam. The response of each pixel to the known number of photons 
constitutes an absolute end-to-end calibration for those pixels. A calibrated 
laser probe is used to measure the number of photons in typical laser pulses in 
the field calibration. A nitrogen laser (337 nm) at a distance of 4 km is well 
suited for such purposes, because the scattered light flux is in the correct range 
and uncertainties due to aerosols are minimal. The scattering angles are just 
greater than 90°, where the differential scattering cross section is minimized 
for aerosols. Moreover, at a distance of 4 km, the extra scattering from the 
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beam by aerosols is approximately canceled by the extra aerosol attenuation of 
the primarily Rayleigh-scattered light flux. Therefore, the roving laser provides 
an independent and redundant absolute calibration of the FD cameras, based 
on the measured light injected into the atmosphere. Using the recorded laser- 
induced ADC traces and the drum calibration constants, we can then compare 
the number of predicted to measured photons. 




Figure 23: Roving vertical nitrogen laser hardware and control system. 

Schematics of the hardware and controlling electronics for the roving laser 
system are shown in Fig. 1231 The roving laser system presently in use is based 
on a 337 nm nitrogen laser, providing about 100 /iJ per pulse, while the output 
is inherently unpolarized. The RMS values of 100 shots in a set, fired over 200 
seconds, is typically 3%. Variations of average energy in a 100 shot group is 
2% on the time scale of 15 minutes, and 10% from night to night, as measured 
with the energy probe inserted into the beam for 100 shots before and after 
the calibration shots. A self-leveling He-Nc laser provides a reference beam 
within 0.03° of vertical, allowing alignment of the nitrogen laser to within 0.1° 
of vertical. 

An embedded Linux PC, communicating via serial connection to a pro- 
grammable GPS timing module, controls the laser firing time relative to the GPS 
second. The time stamp is used to identify calibration events for reconstruction. 

For comparison with the 375 nm drum calibration, the measured relative 
FD response at the two calibration wavelengths must be known. The relative 
response has been measured, and is described in the next section. 

Overall uncertainties of the roving laser calibration have been assessed at 
12%, dominated by laser probe calibration and atmospheric effects. 
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Figure 24: The results of the multi-wavelength measurements (see text, and |37|| .) 

5.6. Multi-wavelength calibration 

For calibration at wavelengths spanning the FD acceptance, a xenon flasher 
is mounted at the back of the drum, with a filter wheel containing 5 notch 
filters for selection of wavelengths. The xenon flasher [3(| provides 0.4 mJ 
optical output per pulse covering a broad UV spectrum, in a time period of a 
few hundred nanoseconds. A focusing lens at the filter wheel output maximizes 
the intensity through the filter wheel and into the light pipe. 

Relative drum intensity measurements at wavelengths of 320, 337, 355, 380 
and 405 nm have been made with the same reference PMT used in the absolute 
measurements. At each wavelength the recorded response from the reference 
PMT, combined with the PMT quantum efficiency and corrected for light source 
and filter width effects, yield a quantity that is proportional to the number of 
photons emitted from the drum. The FD response detected using the various 
filters, with the drum placed in the aperture, can be combined with the results 



from the laboratory to form the curve of relative camera response |3jj shown in 
figure EU 

The curve in figure [5U is the result of an iterated spline fit beginning with a 
response curve predicted from manufacturer's specifications for FD components. 
The shape of this initial curve is dominated by the FD PMT QE and the UV 
filter transmission. The final fit includes effects of the notch filter transmission 
widths (15 nm FWHM), the reference PMT QE, the xenon light source emission 
spectrum, and the relative drum intensity for each filter, all measured in the 
lab, and the observed FD response to the drum for each filter. The relative 
uncertainty at each wavelength on the curve is 5%. 
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5. 7. FD relative calibration system 

The relative optical calibration system [35| is used to monitor detector 
response and to track absolute calibration between drum calibrations. The 
system is used before and after each night of operation. Three positions (A, B, 
and C) arc illuminated for each camera, monitoring different groups of detector 
components. Light is distributed through optical fibers, from permanently 
installed light sources. All components use quartz optics. 




Source A 
with diffuser 



Figure 25: A schematic showing positions of light sources for three different relative 
calibrations of the telescope. 

The A fiber light source is a 470 nm LED [3^. A 7:1 splitter at the source 
provides light to 6 fibers, running to a 1 mm thick Teflon diffusor located at 
the center of each mirror in the FD building, illuminating the camera face. 
The remaining fiber provides light for an output monitoring photodiode at the 
source. The A fiber LED is driven by a constant current source circuit. Normal 
operation pulses the LED for a series of 60 /is pulses. 

Each of the B and C light sources are xenon flash lamps [3|| . Each source is 
mounted at the focus of a f/1.5 lens, with downstream optics including a beam 
splitter (for a source monitoring fiber), a filter wheel, and an f/2.4 lens focusing 
onto a 7:1 splitter. The 7 fibers run to each fluorescence telescope, and to a 
monitor for each output, as for the A fiber. The B source fibers are split near 
each camera and terminate at 1 mm thick Teflon diffusors located at the sides 
of the camera, with the light directed at the mirror. The C source fibers are 
also split, and terminate outside the aperture with the light directed outwards. 
Tyvek sheets are mounted on the inside of the aperture shutters. The sheets 
are positioned such that they are opposite the fiber ends when the shutters are 
closed, and the diffuse light scattered off the Tyvek enters the aperture. 

The B source includes a Johnson-U filter, approximating the full wavelength 
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acceptance of the fluorescence telescopes. The C source filter wheel containing 
interference filters is centered at wavelengths of 330, 350, 370, 390 and 
410 nm, for monitoring detector stability at wavelengths spanning the spectral 
acceptance. 

6. Performance, operation and monitoring of the detector 

All four fluorescence sites have been completed and are in operation. Los 
Leones has been in full operation since March 2004 and Coihueco since July 
2004. Los Morados began data acquisition in April 2005, and the fourth site at 
Loma Amarilla started its operation in February 2007. 

6.1. Uptime Fraction 

The operation of the fluorescence detector can be characterized by the 
uptime, or the fraction of the total time in which the FD was acquiring data. 
Anything that disables the measurement is considered as dead time. The main 
contributions to the dead time are the presence of the sun and nearly full moon 
on the sky, poor weather and the presence of the moon in any phase within 5° in 
the FOV of a telescope. The position of the moon can be calculated in advance 
and the shutters of individual telescopes are closed when the moon approaches. 

The shutters are also automatically closed when the weather conditions 
become dangerous for operation (high wind speed, rain, snow, etc.) and when 
the observed sky brightness (caused mainly by scattered moonlight) is too high. 
The influence of weather effects depends on the season, with the worst conditions 
typically occurring during Argentinian summer. Other periods of dead time are 
caused by the activity of atmosphere monitoring instruments, mostly by lidar 
stations and the Central Laser Facility, readout of the electronics, and any 
hardware or software problems. 

The value of uptime fraction has been derived from the data and is cross- 
checked by several techniques. The average uptime fraction for the whole 
observation period is around 13 % of the total time since operations commenced. 
Averaged uptime fractions for individual telescopes arc shwon in Fig. 1261 For 
Los Leones and Coihueco the values refer to the period January 2005 to January 
2008, for Los Morados and Loma Amarilla the average is calculated from the 
individual start of operation till January 2008. 

The Loma Amarilla building does not yet have a dedicated power line. The 
site is powered by a generator, which is less reliable and has caused a lower 
uptime compared to other FD sites. 

6.2. Background conditions 

The presence of the moon above the horizon increases the background light 
level, which has to be monitored. The direct current induced by the background 
is eliminated by the AC coupling of the PMT base, but it is possible to determine 
background levels using the direct relation between the fluctuations of the sky 
background and the photon flux. The analysis of the fluctuations in the ADC 
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Figure 26: Uptime fraction for fluorescence telescopes. Telescopes are numbered as follows: 
Los Leones (south) site 1-6, Los Morados (east) site 7-12, Loma Amarilla (north) site 13-18, 
Coihueco (west) site 19-24. 



signal (variance analysis) performed for each night of data taking is used to 
monitor the FD background signal and data taking conditions. Thus the amount 
of light is derived from ADC variance value and it is given in units (ADC 
counts) 2 . 

The total background signal is the sum of the electronics background (pho- 
tomultiplier and electronics noise) and the sky brightness (airglow, moonlight, 
stars and planet light, zodiacal light, twilight and artificial light) (2?]]. The 
typical values of background signals arc: 3-5 (ADC counts) 2 for electronics 
background, around 20 (ADC counts) 2 for cloudy nights, and between 25 and 
60 (ADC counts) 2 for clear moonless nights. During nights when the moon is 
above the horizon, the light background level can reach several hundred (ADC 
counts) 2 . 

The optimal background conditions for observation range from 25 to 60 
(ADC counts) 2 which corresponds to photon background flux from approxi- 
mately 100 to 250 photons m~ 2 deg -2 /zs -1 . Under such conditions, one event 
per 2 hours is recorded for which the energy is determined to within 20 % and 
the depth of maximum measured with an uncertainty of better than 40 g/cm 2 . 
These high-quality events are the ones used for physics analysis. 

6.3. Standard Operation 

The FD operation is not fully automated and at present the assistance of a 
shift crew of at least two people per night is necessary. Their responsibilities 
consist of several activities before, during and after measurement each night. 
These include relative calibration of the cameras and optical components 
before and after observation, starting and stopping data taking according to 
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weather conditions, prompt correction of software or hardware defects, etc. The 
operation of the FD is still evolving, and the software development is steadily 
transferring the responsibility of the human crew to automatic operation. The 
ultimate aim is to operate FD telescopes in fully automatic and remotely 
controlled mode from abroad. 

FD is operated in nights with moon fraction below 60 % beginning at the 
end of the astronomical twilight till the beginning of the next astronomical 
twilight. These criteria have evolved over time; before January 2005 the 
maximal illuminated moon fraction was 50%, and hence data taking period 
was 2 days shorter per month. The observation period lasts 16 days per month, 
with an average observational time of about 10 hours (ranging from about 14 
hours in June to 5 hours in December). 

6.4- Slow Control System 

The fluorescence detectors are operated from the central campus in Malargiie, 
and are not operated directly from the FD buildings. Therefore, the main task 
of the Slow Control System (SCS) is to ensure a secure remote operation of the 
FD system. The SCS works autonomously and continuously monitors detector 
and weather conditions. Commands from the remote operator are accepted 
only if they do not violate safety rules that depend on the actual experimental 
conditions: high-voltage, wind speed, rain, light levels inside/outside the 
buildings, etc. In case of external problems such as power failures or 
communication breakdowns the SCS performs an orderly shutdown, and also 
a subsequent startup of the fluorescence detector system if the conditions have 
changed. If parts of the SCS itself fail, the system automatically reverts to a 
secure mode as all potentially critical system states (open shutters, high- voltage 
on, etc.) have to be actively maintained. 

To ensure reliable supervision and to allow for high flexibility and stability, 
the SCS is based on industrial PROFIBUS components. This bus system 
consists of several bus terminals with specific functions, such as analog input, 
digital output, relays, etc. The terminals make up a modular bus-system which 
is addressed and controlled from a PC. The slow control PC is the central 
instance of the SCS in each FD building and runs the main control software 
under a Windows operating system. The system for one of the fluorescence 
detector buildings is sketched in Fig. [23 

The shift crew interacts with the control-system via web-browsers. One 
central webserver communicates via OPC-XML gateways with 4Control OPC- 
servers on the control-PCs. Several views with different levels of detail 
display the status, from an overview of the whole fluorescence detector of the 
Observatory down to single telescopes (see Fig. |2"5)) . 

For the further automation of the data taking, an interface between the Slow 
Control and the Data Acquisition systems is under development. The final goal 
is an automatic, scheduled operation where the shifters have to react only to 
malfunctions of the system. 
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Figure 27: Schematic view of the slow-control system 
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Figure 28: SCS display of a FD telescope 



6.5. Monitoring 

FD data-taking can only take place under specific environmental conditions. 
The light-sensitive cameras must be operated on dark nights with low wind and 
without rain. This makes the operation a full task for the shifters, who judge 
the suitability of operations on the basis of the information given by the SCS. 
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The telescope performance must be monitored constantly to assure the quality 
of the recorded data, as well as guarantee the safe operation of all detector 
components. A user-friendly monitoring tool has been developed to support the 
shifters in judging and supervising the status of the detector components, the 
electronics, and the data-acquisition. 

The monitoring tool [IH is built upon online MySQL databases, organized 
by FD building. The databases record data from regular FD measurements 
and from calibrations and atmospheric surveys. The data are transferred to 
a central server at the Malargiie campus using the internal MySQL database 
replication mechanism. This mechanism recognizes communication problems 
and tries to synchronize database changes when the connection is reestablished; 
this guarantees the completeness of the dataset on the central server, even if 
the information are not available online immediately due to network failures. 
Fig. [29] shows the schematic layout of the databases. 

The user interface is based on a webserver running Apache. The website uses 
PHP, CSS and JavaScript; graphs and custom visualizations arc implemented 
using the JPGraph package for direct PHP calls to the database. They are 
dynamically generated to be accessible not only for the shifter, but also remotely 
for experts. Graph explanations as well as troubleshooting tips are available 
through linked wiki pages. Alarms are triggered automatically in case of special 
occurrences. This tool helps the shifter to monitor constantly the performance of 
the detector to assure the quality of the recorded data as well as guaranteeing 
the safe operation of all detector components. In addition it offers a unique 
opportunity to monitor the long term stability of some key quantities and the 
data quality. 

6.6. Online events 

Besides using the online monitoring system for a continuous performance 
check by integrating over the most recent history of data, the operators are 
furnished with an online event display of the cameras. Two typical event 
examples corresponding to a real cosmic ray air shower and to a background 




Figure 29: Organisation of the databases: The single databases at each FD-building are 
replicated to the database server at the central campus, while other sources like the SD insert 
data directly into the central database. 
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Figure 30: Screenshot of the FD monitoring web interface, showing a selection of FD 
background data for all pixels in the six cameras at Loma Amarilla. 



event are shown in Fig.[3TJ The light spot of the air shower, shown in the upper 
part of the figure, takes about 10 /us to cross the FOV of the camera while 
individual pixels are illuminated for up to about 1 /lis. This is clearly different 
from the event shown in the bottom part of the figure. Here, all activated pixels 
show huge narrow pulses all starting in the same sample of the 10 MHz ADC. 
Such a geometrical pixel pattern and timing feature is incompatible with an air 
shower and is most likely due to a muon penetrating the camera. 

7. Event Reconstruction 

7.1. Geometrical Reconstruction 

A hybrid detector achieves the best geometrical accuracy by using timing 
information from all the detector components, both FD pixels and SD stations. 
Each element records a pulse of light from which it is possible to determine the 
time of the pulse and its uncertainty. Each trial geometry for the shower axis 
yields a prediction for the signal arrival times at each detector component. 

Differences between actual and predicted times are weighted using their 
corresponding uncertainties, squared, and summed to construct a \ 2 value. The 
hypothesis with the minimum value of \ 2 is the reconstructed shower axis. 

In the FD, cosmic ray showers are detected as a sequence of triggered pixels 
in the camera. An example of an event propagating through two adjacent 
FD telescopes is presented in Fig. [32] The first step in the analysis is the 
determination of the shower-detector plane (SDP). The SDP is the plane that 
includes the location of the eye and the line of the shower axis. (See the 
sketch in Fig. 1331 ) Experimentally, it is the plane through the eye which 
most nearly contains the pointing directions of the FD pixels centered on the 
shower axis. (Sec fitted line in Fig. 1321 ) Using a known axis provided from 
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Figure 31: Top: A cosmic ray shower event as it appears in the event display. The pattern of 
the activated pixels is shown in the left panel while the right panel exhibits the response of the 
selected pixels as a function of time for the pixels marked by a black dot. The bin size is 100 
ns. The development of the shower in the atmosphere can be qualitatively seen. Bottom: A 
background event most likely due to cosmic ray muon interacting with the glass of the PMT. 
All activated pixels give signals at the same time, a feature which is not compatible with a 
cosmic ray shower. 



the Central Laser Facility (CLF), described in Ref. the SDP reconstruction 
error can be evaluated by comparing the space angle between the normal vector 
to the experimentally determined SDP and the known true normal vector. This 
uncertainty in the SDP is of the order of a few tenths of a degree depending on, 
for example, the length of the observed track in the camera. 

Next, the timing information of the pixels is used for reconstructing the 
shower axis within the SDP. As illustrated in Fig. [33l the shower axis can 
be characterized by two parameters: the perpendicular distance R p from the 
camera to the track, and the angle xo that the track makes with the horizontal 
line in the SDP. Each pixel which observes the track has a pointing direction 
which makes an angle Xi with the horizontal line. Let to be the time when the 
shower front on the axis passes the point of closest approach R p to the camera. 
The light arrives at the i th pixel at the time 

R 

U = to + — tan[( Xo -Xi)/2]. (1) 
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Figure 32: Light track of a hybrid event as seen by the fluorescence telescopes. The different 
colors indicate the timing sequence of the triggered pixels. The full line is the fitted shower- 
detector plane. (See text for explanation.) The red squares in the bottom left represent the 
surface stations that also triggered in this event. The crosses mark camera pixels that had a 
signal within the time of the trigger, but were marked by the reconstruction algorithm as too 
far either in distance (to the showcr-dctcctor plane) or in time (to the time fit). 




Figure 33: Illustration of the geometrical shower reconstruction from the observables of the 
fluorescence detector [43l |. 



The shower parameters arc determined by fitting the data points to this 
functional form. Using the fast sampling electronics, such a monocular 
reconstruction may achieve excellent accuracy. However, the accuracy of the 
monocular reconstruction is limited when the measured angular speed dx/dt 
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does not change much over the observed track length. An example is shown in 
Fig. [34j For these events (usually short tracks) there is a small curvature in 
the functional form of Eq. JT]) such that there is a family of possible (R P ,Xo) 
axis solutions. R p and xo ar e tightly correlated, but neither value is well 
constrained. This leads to uncertainty in other shower parameters, including 
the reconstructed shower energy. 
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Figure 34: Functional form that correlates the time of arrival of the light at each pixel with 
the angle between the pointing direction of that particular pixel and the horizontal line within 
the shower-detector plane. FD data (color points) and SD data (squares) are superimposed to 
the monocular (red line) and hybrid (blue line) reconstruction fits. The full square indicates 
the SD station with the highest signal. This is a typical event in which the monocular 
reconstruction does not work well. 



The fit degeneracy can be broken by combining the timing information from 
the SD stations with that of the FD telescopes. This is called the hybrid 
reconstruction. The hybrid solution for the example shown in Fig. [M] is shown 
as a blue line and the uncertainties in the parameters are specified in the legend. 

Since the SD operates with a 100% duty cycle, most of the events observed 
by the FD are in fact hybrid events. There are also cases where the fluorescence 
detector, having a lower energy threshold, promotes a sub-threshold array 
trigger (see section 14.2. 2[) . Surface stations arc then matched by timing and 
location. This is an important capability because these sub-threshold hybrid 
events would not have triggered the array otherwise. In fact, the time of arrival 
at a single station at ground can suffice for the hybrid reconstruction. 

The reconstruction uncertainties arc validated using events with known 
geometries, i.e. light scattered from laser pulses. Since the location of the 
CLF (approximately equidistant from the first three fluorescence sites) and the 
direction of the laser beam are known to an accuracy better than the expected 
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angular resolution of the fluorescence detector, laser shots from the CLF can be 
used to measure the accuracy of the geometrical reconstruction. Furthermore, 
the laser beam is split and part of the laser light is sent through an optical 
fiber to a nearby surface array station. Thus, the axis of the laser light can be 
reconstructed both in monocular mode and in the single-tank hybrid mode. 

Using the timing information from the telescope pixels together with the 
surface stations to reconstruct real air showers, a core location resolution of 
50 m is achieved. The resolution for the arrival direction of cosmic rays is 
0.6° These results for the hybrid accuracy are in good agreement with 



estimations using analytic arguments [441 ], measurements on real data using a 



bootstrap method [451 ] , and previous simulation studies [46[ . 

7.2. Shower Profile and Energy Reconstruction 

Once the geometry of the shower is known, the light collected at the aperture 
as a function of time can be converted to energy deposit at the shower as a 
function of slant depth. For this purpose, the light attenuation from the shower 
to the telescope needs to be estimated and all contributing light sources need 



to be disentangl ed [4711: fluorescence light [48|, |49j, |50j, |5l| , direct and scattered 
Cherenkov light |52l. |53| as well as multiple-scattered light 54 , 5|| . 
An example of the measured light at the telescope aperture and the recon- 
structed light contributions and energy deposit profile is shown in Figs. 1351 
and [551 

The calorimetric energy of a shower is estimated by fitting a Gaisser-Hillas 
function [5q | to the reconstructed energy deposit profile and integrating it. 
Finally, the total energy of the shower is obtained by correcting for the 



'invisible energy' carried away by neutrinos and high energy muons 57|. After 
quality selection, the energy resolution (defined as event-to-event statistical 
uncertainty) of the fluorescence detector is < 10% (H)]. 

7.3. Detector Exposure 

The detection volume of fluorescence detectors varies with energy, as showers 
with higher energies emit more light and can be detected further away from the 
detector. The aperture also depends on environmental factors such as night-sky 
background light and atmospheric conditions. To model these effects, large sets 
of detailed MC simulations are used. The response of the Auger fluorescence 
telescopes has been simulated and the detector aperture has been estimated 
as a function of energy, average atmospheric conditions, and primary cosmic 
ray mass for a fixed configuration: a fully built detector with four fluorescence 
detectors and a 3000 km 2 surface array The aperture increases from 



approximately 900 km 2 sr at 10 17 5 eV to about 7400 km 2 sr at 10 19 cV. 

However, the more important quantity here is the hybrid detector exposure 
which accounts for both the time variability of the atmospheric conditions and 
the detector setup. Very detailed time-dependent detector Monte Carlo simula- 
tions are used to reproduce the actual data taking conditions of all components 
of the Pierre Auger Observatory and to derive the hybrid exposure. To assure 
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Figure 35: Example of a light-at-aperture measurement (dots) and reconstructed light sources 
(hatched areas). 
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Figure 36: Energy deposit profile reconstructed from the light at aperture shown in Fig. I35I 
The line shows a Gaisser-Hillas fit of the profile. The energy reconstruction for this shower 
was 3.0 ±0.2- 10 19 eV. 



good agreement between data and Monte Carlo, extensive comparisons are 
performed at all reconstruction levels. 

To limit the influence of trigger threshold effects, an energy-dependent 
fiducial volume has been defined and strict quality criteria are applied during 
the determination of the exposure [4l(. For example, only data with a successful 
hybrid geometry reconstruction are selected, and the SD station used for the 
reconstruction has to lie within 750 m of the shower axis. This condition ensures 
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that the probability to trigger at least one surface station is almost equal to one, 
and it significantly reduces the influence of the a-priori unknown primary mass 
composition. The resulting hybrid exposure, after all selection cuts accumulated 
during the 3 years of the building phase of the observatory, increases from about 
50 km 2 sr yr at 10 18 eV to about 550 km 2 sr yr at 10 19 5 eV. 

8. Summary 

The Observatory has been in scientific operation since late 2003, and the 
synergy between the surface array and the fluorescence detector has proved to 
be most fruitful. In particular, the calorimetric FD energy measurements have 
provided a calibration for the high-statistics data set obtained with full-time 
operation of the surface array. 

This paper has described the Auger fluorescence detector system - its 
hardware, performance, calibration, and event reconstruction methods. The 
24 wide-field telescopes view the atmosphere above the entire 3000-km 2 surface 
array. The Auger Observatory has opened a new chapter in cosmic ray physics 
because of both the large number of recorded high energy air showers and the 
quality of the hybrid measurements. The growing data set is expected to resolve 
the important questions relating to the highest energy particles of the universe. 
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